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Abstract 



A simple expression for the zeros of Weierstrass' function is given which follows 
from a formula for relativistic orbits. 

Index Terms — Elliptic functions, Weierstrass' function, Relativistic Orbits. 

1 Introduction 

For elliptic functions, Weierstrass' function p{Q) is a very important tool. It appears in 
countless applications in mathematics, physics, and engineering. There is some interest to 
determine its zeros. This problem is considered as being difficult. A formula can be found 
in [l], see also [6]. However, to quote ([6], p. 105): "There seems to be no simple way to 
express the roots" . 

In this contribution a simple expression is given for the roots of p(0) = 0. The formula 
follows from a formula for relativistic orbits (see {21 and [4J ) , which in turn originated from 
research in the field of Cauer filters [3j. We assume that the reader is familiar with elliptic 
function theory. We use the notation employed by Lawden [5j. 



2 Relativistic Orbits 

In Lawden's book ([5j, pp. 126-129), it has been shown, that the generalisation of the 
Kepler ellipse to the relativistic orbit leads to the polar equation 

1 

^ ~ A + B ■ sn^{CQ,k) ' 

where r gives the distance of two bodies and Q the angle in the plane of motion with 
suitable reference. As Lawden was mainly interested in deriving the perihel shift of the 
planet Mercury, he only determined approximations of the constants A, B, and C. 
In [2] and [3] these constants have been determined exactly by explicitly solving the orbit 
differential equation for u = 1/r given by 



u" + u = a + 3(3 ■ v!^ . 



The dash denotes derivation with respect to the angle G. The constants a and /3 depend on 
the motion to be considered (e.g. planetary motion, photonic motion, motion of elementary 
particles, etc.). This led to 



B 



c = • ( r. ■ (1) 



Note that the constant C is invariant to the transformation A;^ — t- /c'^ = 1 — /c^. 
The shape for the relativistic motion (and for the magnitude of Cauerfilters) has two 
relevant constants, the value of k and the value of the product a/3. These two parameters 
essentially determine the trajectories (up to scaling). 




3 Simple Expression for Zeros of p 



The orbit equation can be easily solved for the angle = (1/C) • sn~^{y ^^^^^ )- We now 
show the relation of the orbit equation to the Weierstrass function p. To this end we use 
the relation of p to the sn^-function (see f^, p. 553) 

or 

= ^^^) 3 ' 

i.e. C = \/ei — 63 and 63 = —^-^C'^. The value C is given by equation ([T]). The angle Q 
then can be expressed by the integral 

'P^^) dy 



'00 V^y^ - 92y - 93 

where 52 and 53 are the Weierstrass invariants which follow from the differential equation 

p'2 = 4p3 _ ^2 ^ _ 

and can be expressed as sum over all non-zero reciprocal 4-th and 6-th powers respectively 
of the lattice points of the adjoint parallelograms. The polynomial %3 — g^y — g^^ factors 
as 
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= 4(y - 6i)(y - 62)(y - 63) . 




C^, and 61 -|- 62 + 63 = we get 
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which leads to 



53 = 4 6162 63 gs = ^{k^ + l){k^ -2){2e -1)C^ 

52 = -4(6162 + 6163 + 6263) =^ 52 = ^(^^ - fc^ + 1) C"^ • 

Using equation ([T|) we obtain the remarkably simple result 

52 = ^ - a/3 , (3) 

and 

_ (fc^ + l)(fc^-2)(2fc^-l) / l-12a/3 

To obtain a compact description of all cases we compute the discriminant of the poly- 
nomial 4y3 — g2y — 53 which is given by §2 — 27g^ = 16 (61 — 62)^(61 — 63)^(62 — 63)^, thus 
g3 _ 27gi = 2^(1 - k'^fk^C^^ which leads to 



52 - 2751 



(1 - kY -k^-jl- 12af3)' 
28 . (fc4 _ A;2 + 1)3 



If the discriminant is greater than zero we have three real roots, if it is zero we have three 
real roots which are not all distinct and for negative discriminant we have one real and 
two complex zeros. The absolute invariant 52/(2753) follows as 

5i _ {k'-k' + if 



27 gi (A;2 + 1)2 . (^2 _ 2)2 . (jfe2 _ 1)2 • 

The absolute invariant transforms the three cases of the discriminant to the cases that 
52/(2753) is greater, equal or smaller than 1. The absolute invariant is an extremely useful 
tool for characterizing the ratio of the two periods of the corresponding elliptic function. 
It is immediately seen that the absolute invariant is left unchanged by the transformations 
fc2 ^ = 1 - A;2 (like C) and k"^ ^ l/zfc^. 
The zeros of p now follow from equation ([2]) 



- ^ Qo = ±—-sn U ,n ,fe) . 




sn^iCQ) 3 ^ C '\ 1 + k^' 

hence the zeros of the Weierstrass function are given by the simple expression 




, 4(fc4 - fc2 + 1) 

We give some examples. 

Example 1 Computing the length of the famous Lemniscate, one encounters the integral 
J ^^^^4 which can be transformed to the integral J —j^==. Hence 52 = 4 and 53 = 0. 

The values of k'^ follow from equation The possible values are k"^ = —1, k"^ = 1/2, 

and k^ = 2. We select the value of k"^ in the interval [0, 1], i.e. k'^ = 1/2 and obtain 



The approximate numeric value in the parallelogramm around zero is given by ~ 
±1.3110287771 • (1 - i) Clearly this is the value ■ {1 - i). Note that for k = i 

we have Gauss' Sinus Lemniscatus sl{z) = sn{z,i) = ^ sd{\/2z, -^). 

In table I data for some values of k used in the following examples are collected. All 
numeric values for Qq given in the examples are for the parallelogram around zero. 

Example 2 Let §2 = 7, gs = 3, and k = (see table I). This leads to 

Qo = ±y^ • sn-'^i^, PS ±(1.0496381 - i ■ 0.77781243) . 



Example 3 Let 02 = 11, 03 = 7, and k = -^^^-4 



This leads to 



Qo = ±\/2 \/l7- I2V2 • sn-\J- + V2, ^ — -) ^ ±(0.9270373 ± i ■ 0.6766441) . 

V 2 y 2 ± 1 

Example 4 Let g2 = 15, g^ = \/2 • 7, and k = \pl — 1. This leads to 

Go = ±y^ V^3-2\/2 • sn-^(^\/2±\/2, ^2 - 1) f« ±(0.86473386 ± i ■ 0.637892607) . 
We may also use the property 



to obtain additional solutions from known values of 32 and 53 for any A 7^ 0. 

Example 5 Using A = 2^^^'^, we get the zeros of the Weierstrass function for the case 
52 = 15/2"'^/^ and 53 = 7 from the previous example as 

Go = ±2^/12 . /2 4/3 _ 2^ . ,„-i(^y'2 + ^/2, ^/2 - 1) . 

V O z 

In general for given 52 and 5(3 we can easily compute k from equation Using the 
substitution ^ = fc^ ± l/Zc^ — 1 we get a polynomial of degree three in ^ which can be 
solved for ^ from which k'^ and A; follows. Namely we get the polynomial 

f + a-^-a = 
where « = ^ • t^^tt^ — ? > which solved for ^ yields 



with p G {l,e*2'^/3^gi47r/3|^ 'pj^gj^ values of fc^ follow: 

1 



fc2 = 1 . (^ + 1 ± ^^2 + 2e - 3) . (6) 



^Note that a is simply related to the j-invariant 1728 32/(92 — STgl) 



If a real solution for /c^ in the interval [0, 1] exists it is obtained by selecting p = 1 and the 
minus sign in equation ([6]). 



Table I: Data for some invariants 
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4 Conclusion 

The simple and easy to compute expression 

,/ 4(fc4 -k^ + lj _ I 3 

®°=^v — — " ^ 

for the zeros of the Weierstrass elliptic function p has been derived, where k is given in 
closed form as function of (72 and . 
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